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Abstract. Transition towards spatio-temporal chaos in one-dimensional interfacial patterns often involves 
two degrees of freedom: drift and out-of-phase oscillations of cells, respectively associated to parity breaking 
and vacillating-breathing secondary bifurcations. In this paper, the interaction between these two modes 
is investigated in the case of a single domain propagating along a circular array of liquid jets. As observed 
by Michalland and Rabaud for the printer's instability jj, the velocity Vg of a. constant width domain is 
linked to the angular frequency ui of oscillations and to the spacing between columns Ao by the relationship 
Vg — aXow. We show by a simple geometrical argument that a should be close to 1/n instead of the initial 
value a — 1/2 deduced from their analogy with phonons. This fact is in quantitative agreement with our 
data, with a slight deviation increasing with flow rate. 



PACS. 05.45.-a Nonlinear dynamics and nonlinear 
47. 20. Ma Interfacial instability 

1 Introduction 

Many studies have been devoted to the dynamics of one- 
dimensional patterns. One of the motivations is to seek 
for an equivalent of the transition to turbulence in fluids 
(see 12] for a recent exhaustive report). CouUet and looss 
P] have classified the possible behavior in terms of sec- 
ondary instabilities linked to broken symmetries. Most of 



dynamical systems - 47.20.Lz Secondary instability - 

them have been observed in various experimental systems: 
Rayleigh-Benard convection 4 , directional solidification 
IHEEI, Taylor-Dean flow 0, Taylor-Couette flow di- 
rectional viscous fingering between two eccentric cylinders 
("printer's instability") |1I1(JI1H . These behaviors are also 
recovered in numerical investigations of partial differential 
equations governing interface instabilities ^5] > and also by 
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phenomenological models coupling a base mode k and its 
first harmonic 2k I17|. 

Our group has investigated another experimental sys- 
tem: a one dimensional array of liquid columns formed be- 
low an overflowing circular dish (see fig-H'^')- This system 
exhibits dynamical behavior |1HI14I15] similar to those ob- 
served in directional solidification or directional viscous 
fingering. Typical examples of spatio-temporal diagrams 
are displayed on figs. |21 and 13 As in ref. jl3ll4| . these 
diagrams were obtained from pictures of the dish taken 
from above (see insert of fig.H"^)! s-^d by recording grey 
levels along a circle intersecting every column trace. Time 
runs from top to bottom, space (position along the dish 
perimeter) is plotted on the horizontal axis. 




I m ■ 



ffl 



(b) 



Fig. 1. (a) The circular fountain experiment : a pattern of 

Depending on the selected conditions (flow-rate, ini- liquid columns is formed below an overflowing circular dish, 

tial columns mrniber, possible initial imposed motions). Insert : viewed from above, columns appear as U-shaped spots, 

different regimes occur, the most spectacular one being (b) Like in the printer's instability it is tempting to identify 

the spatio-temporal chaos in fig. [2}b. This behavior in- the column array to a chain of springs and beads, 
volves a complex interaction between two elementary de- 

grees of freedom. The first degree of freedom is called a j^^^ temporal chaos. Therefore a careful examina- 

"vascillating-breathing" mode (fig. Ha): adjacent columns ^ion of this interaction in our weU controled system is an 

oscillate in phase-opposition, which doubles the spatial pe- essential step in the study of the transition towards spatio- 

riod but preserves the (x, -x) symmetry The second de- temporal chaos in an extended geometry On the other 

gree of freedom is characterized by domains of asymmetric ^^^^^ spatially chaotic states are so complex that informa- 

cells that break the {x,-x) symmetry. This induces a drift tj^j^g ^aken from regular behavior are to be prefered at the 



of the pattern |13I14I15| . 



present stage. In various systems, including our fountain 



These two bifurcations have motivated several experi- experiment, the vacillating-breathing mode is known to 
mental |5I6I11I13I14II5| as well as theoretical studies |31 1 21 1 6li gi 'IM lift nv the propagation of a solitary drifting domain. 
Among the latter, Misbah and Valance ^3 noticed that which trailing edge is often followed by transient oscil- 
for small systems, the interaction between these two modes lations. The interaction between a solitary parity broken 
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Fig. 2. Spatio-temporal diagrams (I), (a) An oscillating regime 
(VB) extended on the whole dish (F — 0.31 cm'^/s, duration: 20 
s). (b) Regime of spatio-temporal chaos. Interactions between 
oscillations and drift sustain disorder {F = 0.39 cm^/s). 




(b)*t: 



Fig. 3. Spatio-temporal diagrams (II). (a) and (b) Transient 
oscillations following a dilation wave (PB mode) (duration = 
15 s). At low flow-rate {F = 0.11 cm^/s), the VB mode is 
strongly damped (a). At higher flow-rate(r' = 0.29 cm^/s), 
its lifetime reaches one period of rotation of the domain (b). 
Inserted, a magnified view. The dotted lino linking the maxima 
of column oscillations indicates the deviation to a perfect phase 
opposition between neighbouring columns. 



domain and its own oscillatory wake is thus the subject of 
this article. 
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2 Linking oscillation and drift: the phonon 
analogy 

The idea of a possible interaction between oscillations and 
drift in ID patterns is natural, but remains undiscussed 
at least quantitatively by available theories. For instance, 
models based upon the coupling between a base mode and 
its second harmonic ("k-2k approaches") are able to cap- 
ture both behavior analytically 16 but, to our knowledge, 
a solution involving a wall separating oscillations and drift 
has never been proposed. The other well known approach 
based upon symmetry considerations |3ll8j leads to sets of 
coupled phase and amplitude equations that are specific to 
each state (oscillations and drift) and are therefore unable 
to handle their interaction. L. Gil recently showed that 
solitary drifting domains followed by an oscillatory wake 
can be recovered when spatial phase shifts between a base 
mode and a bifurcated oscillatory mode are considered. 
However, this study remained only qualitative. 

On the other hand, quantitative evidences of such an 
interaction were reported by Michalland and Rabaud for 
the printer's instability and confirmed by further inves- 
tigations on the fountain experiment |13ll5j . These works 
all mention a relationship linking the velocity of the do- 
main walls Vg to the angular frequency lu of the tran- 
sient oscillations left behind and the spacing Aq between 
columns (defined outside of the propagating domain) : 

Vg = aXoto (1) 



of a cellular dissipative structure 

where a was found in a range from 0.36 to 0.4. To in- 
terpret qualitatively this result, Michalland and Rabaud 
PP suggested an analogy with the propagation of dilation 
waves on a lattice of springs and beads (fig. ^b). This 
" phonon approach" may seem particularly relevant in the 
case of liquid columns since this pattern appears as com- 
posed of discrete localized singularities jl5| . The numerical 
simulation of spring and mass lattice submitted to a sud- 
den motion of a boundary indeed reveals the propagation 
of dilation waves that exhibits a structure qualitatively 
similar to that of our drifting domains, although they are 
damped by dispersion on long time scales (fig. El-a and b) . 
The velocity of these " domains" and the angular frequency 
of oscillations are linked by a relationship identical to (QJ , 
where Aq is the spacing between each mass, but with a = 
1/2. 

Though this was not explained in details in 1 , this 
value of 1/2 can be deduced from the dispersion relation 
of phonons uj ~ 2{K/my^^ sin |/ca/2| governing the propa- 
gation of waves on the lattice Xn = expi(A;na — wt) (where 
k is the wavenumber, a = Xq the lattice spacing, K the 
spring stiffness and m the mass of each "atom"). Within 
this framework, the domain velocity Vg is identified to the 
group velocity of phonons defined in the long wavelength 
limit (fc — !■ 0), while the frequency of oscillations is identi- 
fied at the boundary of the first Brillouin zone {k = tt/o), 
a limit which corresponds to spatial-period doubling. 

This approach has the merit to provide a simple in- 
terpretation of (p3l- However it is presumably incorrect, 
since a = 1/2 is inconsistent with the measured values 
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to 0.41. As recognized in |15|. including non-linearities phase velocity of oscillations = t^/fc, where k, the 
and dissipation in an improved "phonon" model does not wavenumber of oscillations, is equal to tt/Aq when neigh- 
solve this problem. This is puzzling because Michaland bouring cells oscillate in phase opposition. A similar idea 
and Rabaud argument is rather natural and very general, is also perhaps underlying phase matching conditions in- 
Qualitatively, we think that a different coefficient is voked to describe defect motions in cellular patterns |20|. 
dictated by a presumably non- linear effect: a phase- matching 

condition between oscillation and drift. 4 Experimental Setup 

The experimental setup is simple : silicon oil (viscosity 
3 Geometrical argument ^ = lOO cP, surface tension 7 = 20 dyn/cm, density 

p = 0.97 g/cm^ at 20°C) overflows from a horizontal cir- 

A close examination of spatio-temporal diagrams (fig. OJ 

cular dish (fig. ^a.) at constant flow-rate per unit length 

shows that there is a strong tendancy for phase oppo- 

r. For r between 0.05 cm^/s and 0.6 cm^/s, the liquid 

sition between nearest neighbours in the transient oscil- 

self-organizes as a pattern of liquid columns. This pat- 

latory wake. The phonon model does not impose such a 

tern results from a Rayleigh- Taylor instability combined 

condition. We think this is the only reason for the discrep- 

with a constant liquid supply. It may bifurcate towards 

ancy and postulate that columns oscillate in perfect phase- 
secondary instabilities, which lead to the phenomenology 

opposition with nearest neighbours in the oscillatory wake 

displayed on figs. (21 and O 

as usually assumed for a pure vascillating-breathing mode 

As explained above, these spatio-temporal diagrams 

|3I12| . An idealized sketch of the column motions involved 

are built from pictures taken from above by a video cam- 

at the rear of a drifting domain is sugested on fig. ^a . 

era, the fountain being lighten by a circular neon tube 

In order to obtain a stationary structure propagating uni- 

(see insert of i^g.^-a))- Columns then appear as U-shaped 

formly with time, half a period of oscillations Tr/w must be 

spots, moving along a circle of radius R ~ 4.8 cm slightly 

equal to the time it takes to the rear wall of the drifting 

smaller than that of the dish (5 cm here). The diagrams 

domain to cover the spacing Aq between two oscillating 

are obtained by recording grey levels along this circle. 

columns Xo/Vg. This geometrical argument leads to: 

Thanks to capillary effects, the columns can be manip- 
ulated. One can adjust their number and their initial mo- 

Vg^-Xooj (2) . . 

Ti" tion just by touching them with needles. The rapid coa- 

and hence to a ~ Our argument can also be stated lescence of several columns with another one moving at 

as follows: since trajectories in the (x, t) plane should be constant speed induces a locally heterogeneous pattern, 

continuous, Vg must be equal (in absolute value) to the in which a "dilation wave" followed by damped oscilla- 
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Fig. 4. Idealized geometry involved at the rear front of a prop- 
agating parity broken domain, where the drift has to match 
with transient oscillations. 

tions (or more rigorously a drifting parity broken domain) 
propagates around the dish. 

As shown on and b, the length of the oscillatory 
wake increases with flow-rate which makes possible ac- 
curate frequency measurements. We have systematically 
studied such domains, varying the flow rate F per unit 
length, and we have investigated the evolution of their 
wall velocity Vg, as well as that of the angular frequency 
uj and wavelength Aq left behind. 

5 Experimental results and discussion 

The velocity of domain walls Vg is plotted versus flow- 
rate on fig. [Sj-a. Data are well fitted by a power law, with 
an exposant close to 0.5, which suggests a relationship : 



Vg ^ vr. 



Measurements of the angular frequency lu of transient 
oscillations left behind the propagating domains are plot- 
ted on fig. IHl-b versus flow-rate per unit length F. At a, 
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Fig. 5. (a) Measurements of domain walls Vg versus flow-rate. 
The dotted line provides a fit with power law, with an exponant 
of close to 0.5. (b) Angular frequency of transient oscillations 
left behind a drifting domain (open symbols) and in an ex- 
tended oscillating state (black symbols) versus flow-rate per 
unit length. The flt traces a square root law with a thresh- 
old, (c) Wavelength Aq outside a propagative domain, versus 
flow-rate, for several domain size. 
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few percent of accuracy, this frequency is very close to uation a = I/tt implied by the argument suggested on fig. 

that of global oscillations (black symbols) such as those EJ-a only occurs close to the limit flow-rate Fc = 0.06±0.01 

on fig. [21-a. As in previous works |13I14I15] . uj increases cm^/s. It is interesting to note that this limit flow rate 

with r . nearly coincides with the critical threshold of appearance 

Figure El-c shows measurements of the wavelength out- of oscfllations (fig.|5}b). Indeed, angular frequency mea- 

side a propagative domain (Ao). This length does not vary surements are correctly fitted by the following law: 
with flow-rate, neither with domain size. Its value is around 



r — r 

1.08 cm. ijj = ijjQ + (4) 



D 

This is consistent with a Hopf bifurcation. If this really 
holds, it would suggest that only holds in the limit of 
a vanishing oscillation amplitude, the deviation from I/tt 
being intrinsically non-zero for any flow rate. 

Finding a = I/tt instead of 1/2 shows that the liquid 



Despite the obvious limitations of such an approach, it 
is here interesting to confront this result to the "phonon 
analogy" . In this point of view, assuming an effective col- 
umn mass proportional to F, this analogy would imply 
that the effective stiffness K could scale as . This depen- 
dance is compatible with an interaction between columns 

column array and an equivalent spring and beads sytem 

through inertial terms of Navier-Stokes equation, the Reynolds 

are rather different in the details although the analogy 

number Re = F /v being here of order 1. 

In fig. Oa, we have plotted the values of Vg versus 
Xqlo. The slope appears to be very close to I/tt, at low 
flow-rates, in agreement with eq. (|2J), and seems to be- 
come slightly larger when the flow-rate increases. In fact 
the situation is more complicated. As appears on flg. El-b, 
there is a systematic difference between Vg and Aqw, that 
increases linearily with _r, following the empirical law : 



can perhaps be useful to get reasonable scaling laws. As we 
have explained, the value found for a seems to be dictated 
by : (1) a strong tendancy for neighboring columns to os- 
cillate in perfect phase opposition and (2) a "coherence" 
condition between the wall motions and the oscillations. 
In terms of Fourier analysis the first condition can be ex- 
pressed by fc = tt/Aq, i.e. the system lies at the boundary 
of the first Brilloiun zone. 

The second condition reads Vg = = u/k, i.e. the 
Xquj tt D ^ ' velocity of the domain walls must be equal to the phase 

D has the dimension of a diffusion coefficient. Its value velocity of the oscillations. The first condition is rather 
is around 4.6 cm^/s. The threshold F^ is around 0.06 reasonable if we think in terms of secondary bifurcations 
cm^/s. of a dissipative pattern. Near the threshold of oscillations. 

The effective value defined by ae// — Kj/Aqw varies we can speculate that a single wavevector and a single 
from 0.31 to 0.38 linearly with flow-rate. The "ideal" sit- frequency are selected. In some sense this mode is the 



1-^-^^ (3) 
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Fig. 6. (a) Velocity of propagative domains versus the quan- 
tity Xquj with Lu = 2tt/T. (b) Difference between Vg/{Xou;) and 
I/tt versus flow- rate per unit length, (c) Evolution in the plane 
(fc, uj) where k is the wavenumber of the representative points 
of the residual oscillations. Flow rate F increases from right to 



sole excited eigenmode of the problem, whereas in a lat- 
tice of springs and beads, any eigenmode involved in the 
dispersion relation can contribute to the oscillatory wake 
following a dilation wave. Since a differs from I/tt with in- 
creasing flow-rate means that at least one of these two con- 
ditions is progressively relaxed. A careful examination of 
the spatio-temporal diagrams has convinced us that only 
the first condition is violated, the second one holding at 
any value of the flow rate. 

This is clearly visible on the insert of fig.O-b. Isophase 
lines linking second neighbours are slightly inclined, while 
the structure of the trajectories remains consistent with 
the qualitative picture of fig. ^a. To be more quantitative, 
we have plotted onOc the path followed by the system in 
the plane (fc, uj) when F is increased. Here we recall that 
to designates the frequency of the transient oscillations, 
while k is their wave vector (which is different from the 
wave vector of the pattern fco), so that locally the column 
position varies as exjp[i(kx — tot)]. The black symbols are 
obtained by direct determinations of k on spatio-temporal 
diagrams across slope measurements of the isophase line 
linking second neighbors (dotted line on fig. ^a.). This 
slope ip reads: 



tan^j — —{k — n/a) 



(5) 



The open symbols are obtained by using the relation- 



left, from r = 0.076 cm^/s for fc=2. 83 cm ^, to r = 0.36 Cm^/S I- , nr 1 • TA Tr rn, r . , < 

I ' I ship /J = w/V^, and assummg — Vg. ihe fact that both 

for fc=2.31 cm~^. Black circles are based upon direct k mea- 



surements from spatio-temporal diagrams, while open symbols 
stand for values deduced from k = i^/Vg. The two points inside 
the circle correspond to negative values of Vg/XoL^i — I/tt. 



kind of symbols overlap shows that for any value of F, the 
relationship — Vg holds. Since the fiowrate 7^ increases 
from right to left of this graph, the wave number k of 
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perturbations of columns positions is equal to = 2.9 
cm^^ at low flow rate and progressively decreases when 
r increases. In other words, the system initialy lies at the 
boundary of the Brillouin zone, and as F increases, k be- 
comes smaller than tt/Ao, the perfect phase opposition is 
progressively lost. 

6 Conclusions - Conjectures 

In summary, this paper reports accurate measurements 
which allow us to evidence a fundamental relationship 
linking the velocity of a propagating parity broken do- 
main to the frequency of its oscillatory wake. At a few 
percent of accuracy , this pulsation coincides with that of 
the oscillatory mode itself observed alone in an extended 
geometry (fig. El-b), which suggests that it is finaly this 
oscillatory mode that rules the propagation of parity bro- 
ken domain. Similarities and differences with phonons on 
a lattice of springs and beads have been discussed in the 
plane (k,uj), the residual oscillations exhibiting a possible 
shift of the wavenumber k with respect to the boundary 
of Brillouin zone of order of 20 % of the maximal k value. 

We believe that the problem adressed in this article is 
important for several reasons. First, models based upon 
symmetry arguments 3 miss the relationship reproduced 
in eq. In this approach, uj is just a free parameter 
that is selected at will to rebuild the spatio-temporal di- 
agrams starting from the amplitude evolutions. This sug- 
gests that improvements of this approach must be built 
|19|. In another direction, it is to note that k-2k mod- 
els |l(jll7) are able to capture both oscillations and drift. 
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Therefore, a possible promising other idea to interpret our 
result would consist in building a model of the wall sep- 
arating oscillations and drift by starting from Caroli et 
al equations suggested in ref.^T. This could constitute a 
better framework to recover eq^by a rigorous calculation. 

Next, the point investigated here has something to do 
with an insufficiently studied problem, i.e. that of Galilean 
invariance in pattern dynamics. For instance, CouUet and 
Fauve 22 studied the effect of this invariance on a Ginzburg- 
Landau-like equation and discussed the consequences for 
systems in which this invariance is slightly broken by rigid 
boundary conditions. Another paper from Shraiman |21| 
is also available in the case of the Kuramoto-Shivashinsky 
equation. Both papers show that, in such systems, a phase 
dynamics of second order in time should be observed. This 
second order time dynamics is an alternative to the idea to 
invoke a column or cell inertia. Finaly, this subtle interac- 
tion between oscillations and drift is certainly important 
in the genesis of spatio-temporal chaos, because it can in- 
fluence the nucleation process of defects, via "shock" for- 
mations in the phase field. We hope that our paper will 
motivate further studies in this field. 
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